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A formula is derived that allows the computation of one-loop mass shifts for self-dual 
semilocal topological solitons. These extended objects, which in three spatial dimensions 
are called semi- local strings, arise in a generalized Abelian Higgs model with a doublet 
of complex Higgs fields. Having a mixture of global, SU(2), and local (gauge), £7(1), 
symmetries, this weird system may seem bizarre, but it is in fact the bosonic sector of 
electro- weak theory when the weak mixing angle is \. The procedure for computing the 
semi-classical mass shifts is based on canonical quantization and heat kernel/zeta function 
regularization methods. 
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1 Introduction 

The purpose of this paper is to address the computation and analysis of one-loop shifts to the 
masses of semi-local planar topological solitons arising in a natural generalization of the Abelian 
Higgs model. Seen in (3+l)-dimensional space-time these solitons become semi- local strings 
whereas their masses give the string tensions, see PQ-[2]. At the critical point that marks the 
phase transition between Type I and Type II superconductivity, the set of semi-local self-dual 
topological solitons is an interesting 4/-dimensional moduli space, where I is the number of 
quanta of the magnetic flux, see [3]- [I]. Recently, superconducting semilocal (non self-dual) 
strings with very intriguing properties have been discovered in this model [3]. 

Computations of one- loop mass corrections will be performed using the heat kernel/zeta 
function regularization method. The high-temperature asymptotic expansion of the heat func- 
tion, see [Tj-[6]-[H]-[9], is a powerful tool that was applied for the first time to the calculation of 
kink mass shifts in pU] -the M = 1 SUSY case- and pj] -the non SUSY case-. One- loop cor- 
rections to Af = 2 supersymmetric self-dual Nielsen-Olesen vortices were computed in a similar 
approach by Vassilevich and Rebhan-van Nieuwenhuizen-Wimmer in References [12] and [T3] . 
In the second paper, the authors also showed that the central charge of the SUSY algebra is 
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classical level. Sometime later, we calculated the one-loop mass shift for N = (non-SUSY) 
self-dual NO vortices carrying a quantum of magnetic flux in Reference |15j . Mass shifts have 
been given for spherically symmetric self-dual vortices (when several solitons of a quantum of 
flux have coinciding centers) in [16] up to four magnetic flux quanta. Cruder approximations 
were also provided for the mass shift of two separated self-dual NO vortices -each of them with 
a quantum of magnetic flux- as a function of the inter-center distance. 

In Reference [T7] , we studied the one-loop correction to the energy of a degenerate manifold 
of kinks that arise in a very interesting family of models with two real scalar fields in (1+1)- 
dimensions. These field theoretical systems are obtained through dimensional reduction -plus 
a reality condition- of an M = 1 supersymmetric Wess-Zumino model with two chiral super- 
fields, see [H]-|19j-[20j. A comparison between the mass shifts of these composite kinks and 
the correction to the mass of the ordinary A02 kink was offered in [2TJ. In this paper, we 
address a similar, but more difficult, situation in (2+l)-dimensions, comparing one-loop mass 
corrections of the topological solitons that arise in two planar Abelian gauge systems; one with 
two complex scalar fields and the other with a single complex scalar field. The methodology 
used to accomplish this task is explained in detail in Reference [22], where a complete list of 
References can be found. 

Our paper is organized as follows: in Section §2 we describe the model and develop pertur- 
bation theory around one of the vacua. A one-loop renormalization is also performed. Section 
§3 is devoted to summarizing the structure of the moduli space of self-dual topological soli- 
tons. As a novelty, we also apply a variation of the de Vega-Shaposnik method |23j to find 
numerical solutions for spherically symmetric topological solitons. In Section §4, we explain 
how to obtain one-loop mass shifts in terms of generalized zeta functions of the second-order 
differential operators ruling the small fluctuations of the bosonic and ghost fields, and in §5 
the high-temperature expansion of the heat kernel is used to give the final formula for one-loop 
mass shifts of semi-local self-dual topological solitons after application of Mellin's transforms. 
We present our results in Section §6 by means of Mathematica calculations of the coefficients 
of the asymptotic series giving the heat functions. Finally, in the Appendix we offer several 
Tables where these coefficients are shown. 

2 The planar semilocal Abelian Higgs model 
The Model 

The semi-local Abelian Higgs model pQ describes the minimal coupling between an [/(l)-gauge 
field and a doublet $ of complex scalar fields in a phase where the Higgs mechanism takes 
place. The term semilocal refers to the fact that while the global symmetry of this system is 
SU{2) x U(l) only the U(l) factor is gauged. Defining non-dimensional space-time variables, 
x M — > — x M , and fields, $ — > v&, — > vA^, from the vacuum expectation value of the Higgs 
field v and the ?7(l)-gauge coupling constant, e, the action for the semilocal Abelian Higgs 
model in (2+l)-dimensions reads: 

--F^ + \{D^D^ - ^V<t>K) - 1) 2 1 
4 Z o 

where the covariant derivative is defined as = — zA M )$. The physical spectrum differs 
from that of the standard AHM model in that, along with the massive vector boson and Higgs 
scalar, there is a complex Goldstone field. The parameter k 2 = -4 measures the ratio between 
the square of the masses of the Higgs, M 2 = Xv 2 , and vector particles, m 2 = e 2 v 2 . Here, A is 
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of length x mass. 



Feynman Rules in the Feynman-'t Hooft R gauge 

In order to obtain the Feynman rules, we expand the action around a classical vacuum. We 
shall profit from the SU(2) global invariance and we choose a vacuum <& v with real upper and 
vanishing lower components. The shift of the fields 

_ / l + H(x»)+iG(x»' 

corresponds to the identification of H(x fl ), and f(x tJl ), respectively as Higgs, real, and 

complex Goldstone fields. The choice of the Feynman-'t Hooft R-gauge 

R(A^G) = 8^ + G 

requires a Faddeev-Popov determinant to restore unitarity, which amounts to introducing a 
complex ghost field x- All this together allows us to write the action in the form 



S = V - I iV'.r 



+ -d^G - -G 2 + foH&H - y 2 + d^*d»<p 
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- y#(# 2 + G 2 ) + A^HG - d"GH) + H{A^ - x *x) + iA^d^ip - <pd>*i? 

+ A^\<p\ 2 - ^{H 2 + G 2 ) 2 + \{G 2 + H 2 )A^ - ^|^| 2 (|^| 2 + H 2 + G 2 + 2H) 
o 2 2 

in order to deduce the Feynman rules for the propagators and vertices shown in Tables 1 and 
2. The propagator of the complex Goldstone boson plus two trivalent and four tetravalent 
vertices accounting for the interactions of Goldstone-anti-Goldstone pairs must be added to the 
Feynman rules of the Abelian Higgs model in the Feynman-'t Hooft gauge. Nevertheless, in 
order to guarantee that all physical quantities will be ultraviolet finite there is still a remaining 
piece to be added: the action for the counter-terms. We shall compute this up to one-loop 
order in the sequel. 



One-loop mass renormalization counter-terms 

All ultraviolet divergences come from the integral: 

d 3 k i 



(2tt) 3 k 2 - c 2 + ie 



We also note that even if there are massless particles propagating in 2+1 dimensions, 1(0) is 
infrared convergent. The renormalization of the semilocal AHM requires the computation of 
the following one-loop graphs: 

1. Higgs boson tadpole: 



+ - -XJ + - + -XJ + v ) 

= -2i(k 2 + 1) • 1(1) - in 1 ■ 1(0) + finite part 
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Particle Field Propagator Diagram 
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2. Higgs boson self-energy: 



+ : + i * + -TV+ 



= -2i(n 2 + 1) ■ 1(1) - in 2 ■ 1(0) + finite part 

3. Real Goldstone boson self-energy: 

Q + C + + O + " 

= -2i(n 2 + 1) • 1(1) - in 2 ■ 1(0) + finite part 

4. Complex Goldstone boson self-energy: 

+ + + + -o - 

= -Ai(n 2 + 1) ■ 1(1) - 2in 2 ■ 1(0) + finite part 

5. Vector boson self-energy: The potentially divergent part is 



+ -|- ™{ )VW -)- J -\- ryJXfyi VWV 

-VVAAAA/WV '* *WVYA7wW 

= 2i ■ [1(1) + l(0)]g^ v + finite part 

In (2+l)-dimensions the graphs above are the only ultraviolet divergent diagrams for any 
number of loops in the diagrams, not only in one-loop order, and the theory is super-renormalizable. 
Thus, in a minimal subtraction scheme we get rid off all ultraviolet divergences arising in the 
vacuum sector of the model by adding the counter-terms 

£tt. = \ [2(^ 2 + 1) ■ J (!) + ■ • PiM^iM + ^MSaM - 1] 



C£ t = + 1(0)] ■ A,(x»W(x 
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Vertex Weight Vertex Weight Vertex Weight Vertex Weight 




Table 3: Counter-term vertices 



Diagram Weight 




i[2{k 2 + 1)7(1) + k 2 I(0)} 
i[2(k 2 + 1)7(1) + k 2 I(0)} 
i[2(k 2 + 1)7(1) + k 2 /(0)] 
i[A(k 2 + 1)7(1) + 2k 2 I(0)] 



W^VVV -2*[/(l) + /(0)]p- 



which gives rise to the vertices shown in Table 3. In our renormalization scheme, finite renor- 
malizations have been adjusted in such a way that, on one hand, the divergence due to the 
tadpole graph is exactly canceled when the theory is fine-tuned to the so called self-dual limit 
k, 2 = 1 and, on the other hand, the global SU(2) symmetry remains unbroken up to one-loop 
order. 



3 Semilocal vortices 
Finite energy solutions 

Apart from the ground states built around the S 3 classical homogeneous solutions, the semilocal 
AHM has room for other quantum states arising from the extended classical configurations, 
which are stable owing to topological reasons. These configurations are the topological solitons 
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such that their energy, 

£ = J cPx^FaFa + ) t A$ + y ($ f $ - l) 2 ] , 

is finite. The topological character of these solutions is peculiar in that although the scalar 
vacuum manifold is S* 3 , and thus simply connected, the configuration space C 

C = {$(£) g Maps(M 2 ,C 2 ), G Maps(R 2 , TM 2 )/£($, A { ) < +00} 

is the union of Z disconnected sectors. Finite energy configurations show the asymptotic be- 
havior 

$ f $k = l , D^\ sL = (d^-iA^)\s L =0, 

where is the circle that bounds the plane at infinity. Parametrizing by the polar angle 
6 = arctan^, and, modulo the global S , f/(2)-symmetry, choosing = $ y for 9 = 0, the 

boundary conditions on the covariant derivatives provide a map, — > S\, between the sphere 
at the infinity spatial and the fiber Si at the north pole of the base S 2 of the S 3 vacuum Hopf 
bundle: 

Continuous maps between one-dimensional spheres are classified according to the first homotopy 
group and, because the temporal evolution is continuous, n (C) = n^S^) = Z, the zero 
homotopy group of C is non-trivial. Thus, C = U/gzC/ is the union of disconnected sectors 
characterized by an integer number I. 

Moreover, the boundary condition for the vector field A^s^ = —i&di&ls^ is tantamount 

to: 

A t (x) \ SL = ^^(x) \ SL - 2 (f)|^-(f) \ SL + 3 (f)|^i(f) |^ - 4 (f)|^(f) |^ 

/ piie \ 

The 5C/(2)-orbit of = f J is: 

v _ ( e iei sin^ e* e2 cos^ \ ( e iW \ _ ( e'^+^sin^ \ 
1 ~ \ -e-^cos^ e-' l9l sm^ J \ J ~ \ -e^-^cosip J 

where 61,62 € [0, 2n] and ip G [0, |] are the Hopf coordinates of the S* 3 sphere. Therefore, 

36 

A t (x) \ sio = 1-k- (cos 2 (#! + I6)sm 2 ^ + sin 2 (^ + Z^sin 2 ^) + 

OXi 

36 

+ I— (cos 2 (# 2 - 16)cos 2 ^j + sin 2 (# 2 - W)cos 2 ^) 

OXi 

I 89 

dxi 

and the topological (winding) number I has a direct physical interpretation in terms of the 
magnetic flux carried by the planar soliton: 

g = j (A^dx 1 + A 2 {x)dx 2 ) = I j> (j^dx 1 + T^rdx 2 ^J = I jT dO = 2nl 
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We shall restrict ourselves to the critical point between Type I and Type II superconductivity: 
k, 2 = 1. The energy can be arranged in a Bogomolny splitting: 

E = f ^ ± zD 2 ^(D^ ± i£> 2 $) + [F 12 ± i($t$ - l)] 2 ) + 

One immediately realizes that the solutions of the first-order equations 

D 1 $±iD 2 <$> = , Fi 2 ± -($ f $ - 1) = (1) 

are absolute minima of the energy, and are hence stable, in each topological sector that has 
a classical mass proportional to the magnetic flux. Because the first-order equations can be 
obtained from the self-duality equations of Euclidean 4D gauge theory through dimensional 
reduction, the vortex solutions of ([T]) are called self-dual at the limit k 2 = 1. 

We follow [3] to summarize the properties and existence of the so-called semi-local self- 
dual topological solitons: the solutions of (00) for non-negative / (plus sign in the first-order 
equations )0. The equation on the left in (DO) is tantamount to: 

A(z,z) = -id s log ,0 = 1,2 , ) > ( 2 ) 

where the complex notation for coordinates and fields is: 

z = xx+ix 2 , d z = -(<9i - id 2 ) , A = -(Ai - iA 2 ) 
From (j2J), one sees that: 



$1(2,2); v y $i(z,z) Pi(z) 

where Q m and Pi are polynomials of respective degree m and Z in z, such that a; (2) is locally 
analytic. The behavior of the Higgs field at infinity (up to global SU{2) transformations) 
compatible with finite energy requires that m < I and Pi{z) be monic: 

Pi(z) = z l + pi^z 1 - 1 + ■ ■ ■ + Pl z + po , Qi-i(z) = qi-iZ 1 - 1 + ■ ■ ■ + q lZ + q 

Therefore, the moduli space of semi-local topological solitons depends on 21 complex (4/ real) 
parameters: (p a , q a ), a = 0, 1, • • • , I — 1. Varying the values of (p a , q a ) one varies the / zeroes of 
$1 and the I — 1 zeroes of $2 i n such a way that the locations of the zeroes of the two Higgs 
fields plus the scale and orientation of $ 2 parametrize the SSTS moduli space. Equation (pQ) 
on the right becomes 

A u(z, z) + l- e u ^ = A log (\P\ 2 (z) + |Q W | 2 (z)) , (3) 

where 

u{z,z) = n(z,z) + log(l + \uj\ 2 (z)) , tt(z,z) = log ^!(z,z) 



x It is trivial to obtain the solutions for negative I if the solutions for positive I are known; simply, complex 
conjugation gives the solution of equation (TT]) with the other signs. 
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a solution of © compatible with finite energy boundary conditions [Sill]- Therefore, the finite 
energy solutions of ([I]) with a magnetic flux 2irl take the form 



where -u(;z, z) is a solution of ([3]). In the I = 1 case we have that: 



■ e 



e 2 



u(\z-z \; go) 



Z0\ + W\ 

Setting the parameter g to zero, we find the embedded Nielsen- Olesen vortex centered at z : 

$i 1] (z,z;0) \ = f jfz^f ^ _ e |u(|*-«o|; o) 



Alog[|2;-zo| 2 + M 



4| 


go! 


2 


(I* - ^0 


I 2 + 


gol 


2^2 



however, tends to zero for very large |go| and the solution of becomes: u(\z — zq\; |go| » 
1) ~ 0. Therefore, 

$ ( {= 1) (z,z;\q \»l)\ = 1 fz-zo 

is precisely the field profile of the lump centered at zq with radius |go| and topological charge 1 
in the planar CP 1 model. The Higgs fields spreads over the vacuum manifold S 3 for very large 
| go I whereas the ANO profiles are found for very small |go|. Self-dual semi-local topological 
solitons interpolate between self-dual ANO vortices and CP 1 -lumps when |g | varies between 
and oo. 



Self-dual semi-local topological solitons with spherical symmetry 

Among the topological solitons, the simplest ones are those in which all the roots of both 
Pl(z) — (z — zo) 1 and Q m (z) = q m (z — z ) m are at the same point. Let us focus on this case 
by considering the spherically symmetric ansatz around this point chosen at the origin. If 
r = +\/xf + irf, for configurations of the form 

= ( m\k~ » J • /(r) = ^ + ' es ' ' 
MxuX2) = _ lElj ^> Xj , , i(r) = _fc_. e H-..) , 

the first-order PDE's reduce to a system of nonlinear ordinary differential equations: 

\tr = ~^(r) + \Kr)\ 2 -l) (4) 

| = l -f(r)[l-a(r)} (5) 

* = l -\h\(r)(^ - a(r)) (6) 
8 



lim f(r) 

r—*oo 

/(O) = o 



lim h(r) = 

r— >oo 

I MO) I = \ho\S m ,o 



lim a(r) - 

r^oo 

a(0) = 



(7) 



required by energy finiteness plus regularity at the center of the vortex. Let us recall that the 
boundary conditions at infinity also require that mbe < I. The magnetic field and the energy 
density of the spherically symmetric vortices in terms of the field profiles f(r), a(r) are: 

I da 

B[r) = — 

v ' 2r dr 

S(r) = \(y 2 + 1)(1 - f(r) - \h(r)\ 2 ) 2 + - «(r)) 2 + - a (r)) 2 . 



Semi-local strings with a quantum of magnetic flux 

We now go on to the most elementary solutions that carry a quantum of magnetic flux, or, 
I — 1 = to + 1 . We follow the procedure developed in |23j to solve the non-linear ODE 
system fll])-©-© with boundary conditions ((7|). First, we consider small values of r and in 
the first-order differential equations we test the power series 



f(r) = fi-r + f 2 -r 2 + f 3 -r 3 + f 4 r 4 + 



air 



ot\ ■ r + «2 • r 2 + «3 • r 3 + «4 • r 4 + 



h(r) = h + hi ■ r + h 2 ■ r 2 + h 3 ■ r 3 + h 4 ■ r 4 + 



(9) 
(10) 



where /, and ot,-, j = 1, 2, 3, ■ • •, are real, whereas hj, j = 0, 1, 2, • • -, are complex coefficients. 



The coupled first-order ODE's are solved at this limit by (l8l)-(l9l-( jT0l) if 

f(r) 



fi-r+ ^(\h \ 2 - 1) • r 3 + A [(fa? _ i)( 2 |/ i0 | 2 - 1) + Afl] ■ r 5 + . . . 



a(r) 



\(l - \h \ 2 ) ■ r 2 
\h \\\h \ 2 - 



^\ho\ 2 (\ho\ 2 -l) + lf? 



■r 4 - 



768 



l)(Wo\ 2 - 2) + A/i(5|M 2 -4) 



■r° + . 



h{r) ~ 
We stress that 



ho + ^(\h \ 2 - 1) ■ r 2 + A [(i^i- _ 1)(2N 2 _ 1} + 4/ 2 ] •/ + ... 



/in 



9o 
l?o| 



e 2 



5«(0; go) 



is determined by the behavior of the solution at the origin. After that, only a free parameter, 
fi, is left in the exact solution near the origin. Second, a numerical scheme is implemented 
by setting a boundary condition at a non-singular point of the ODE system, which is obtained 
from the power series for a small value of r (r = 0.001 in our case). This scheme prompts a 
shooting procedure by varying f\, where the correct asymptotic behavior of the solutions is 
obtained by setting an optimal value for fi for a given value of h . Finally, the first-order ODE 
system is solved for large r by means of a power series in -: 



3=0 



a[r) 



E 

3=0 



a J ■ r 3 



h(r) 
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f(r) 



air) 



\h(r)\ 



5 

16 



^| 4 )-r- 6 + 



1|2 



■1152 + 158|/i 



1 [2 



r~ s + . . . 



4 1 1 128 1 1 
fe 1 ! 2 ■ r~ 2 + l/iTf-8 + l/i 1 ! 2 ) • r~ 4 + |/i 1 | 2 (-192 + 241/i 1 ! 2 



|^| 4 )-r- 6 



+ |/i 1 | 2 (9216 + 1120I/1 1 ! 2 - 48I/1 1 ! 4 + (/i 1 ! 6 

/i| Q 1 * 1 I 3 / ™ ^ 



r -8 + 



~ \h\ 1 -r~ 1 + 
+ l^l 3 



r~ 3 + \h 
5 



2 + -l/i 1 ! 2 ) -r" 5 + 



32 + 5I/1 1 ! 2 - — l/i 1 ! 4 ) 
ii lg i i ; 



r ' + . . . 



Again, only one free parameter, h 1 , is left. The value of h 1 is fixed by demanding continuity of 
the solution at intermediate distances (r = 15 in our case) obtained by gluing the short-r and 
large-r approximations. In particular, this has the important implication that 



\ho\ = 



\h}\ = 



linking the null value of \h \, which gives the embedded ANO vortex, with the null value of the 
constant \h}\ setting the behavior of the solution for very large r. Another important remark 
is that the large r behavior of self-dual semi-local defects differs from the large r behavior of 
self-dual ANO vortices that decay exponentially. 

The following Figures show the results obtained with this procedure for several values of 
h . Note that ho = ijjr • e^°' qo \ u(0; 0) = -oo, u(0; oo) = 0, ho = for the ANO vortices, 
and ho — 1 for the CP 1 -lumps. It is observed in the graphics that the field profiles reach their 
vacuum values at distances of the order of r = 15. Consequently, almost identical numerical 
solutions would be generated by sewing the numerical and the asymptotic solutions together 
at r greater than 15. 




10 15 20 





a) Functions f(r), h(r) and a(r) and b), c) Energy density for ho = 0.1 






a) Functions f(r), h(r) and a(r) and b), c) Energy density for ho = 0.3 



10 



a) Functions fir), h(r) and a(r) and b), c) Energy density for ho = 0.6 




a) Functions f(r), h(r) and a(r) and b), c) Energy density for ho = 0.9 



Figure 1: Field Profiles and Energy Densities for Semi- local Topological Defects 

4 One-loop correction to the masses of 
semi-local self-dual topological solitons 

Casimir energies 

The starting point for computing the quantum corrections to the soliton mass is the evaluation 
of the zero-point oscillations of the fields around both the soliton and the vacuum. The mea- 
surement of the zero-point soliton energy with respect to the vacuum energy shows a strong 
analogy with the Casimir effect, where the plates are substituted by the soliton profile. 
Let us first consider small fluctuations around the semi-local soliton by writing: 

= S(x) + 5S(x ,x) , Aj(x) = Vj(x) + 5aj(x ,x) 

By S(x) and Vk(x) we respectively denote the scalar and vector fields of the topological soliton 
solution, whereas 5S(xq,x) and 6a,j(xo,x) calibrate the small deviations of the bosonic fields 
with respect to the classical solution. Of course, not all fluctuations are physically relevant, and 
we should avoid pure gauge deformations. In order to do so, we impose the Weyl/background 
gauge condition: 

A (x ,x) = , dj8a,j(x ,x) + -(S^ (x)5S(x , x) — 8S^(x ,x)S(x)) = 

Under these constraints, the ground state energy in the topological sectors up to 0(5 3 ) order 
(one-loop) is: 



ii 



5^(x ,x) 



( 5a 1 (x ,x) 
Sa 2 (x 0l x) 
5Sl(x ,x) 
SSf(x ,x) 
5Sl(x ,x) 

\ 5S%(x ,x) J 



and 5x{x) corresponds to ghost fluctuations. K and K G are second-order differential operators. 
K is matrix- valued and fairly complicated 
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C 



\ 



V 



and K G is the scalar differential operator: 



G 



K 



In all these formulas, the notation is 



S(x) 



Si(x) 
S 2 (x) 



-d k d k + ISA 2 + \S 2 \ 



S{(x)+iSl(x 
Sl(x) + iSl(x) 



and analogously for 5S. The covariant derivatives are VjSfj 
differential operators in the diagonal of K take the form 



f). Q a _i_ F ah V- Q b 

U J°M + t v 3°Mi 



and the 



A = -dA + lS^ + lS.l 2 , B 
C = -d k d k + ^(\S 1 \ 2 + 3\S 2 \ 2 + 2V k V k 



-d k d k + -(3|5x| 2 + |S 2 | 2 + 2V k V k - 1) 



The energy due to the zero point fluctuations around the vacuum must be subtracted. As 
is fairly obvious, the Hamiltonian for the vacuum fluctuations has exactly the same form as 
( TTTT) for the vacuum values: Si — 1, S 2 = 0, V\ = V 2 = 0. We shall denote by K and Kq the 
operators K and K G in the vacuum background. 

When all the fluctuation modes are unoccupied, we formally obtain the contribution to the 
semi-local topological soliton ground state energy as a difference between two "super-traces" of 
differential operators acting on column vectors of L 2 (R 2 ) functions. The first super-trace comes 
from (the star stresses the fact that zero modes do not contribute to the trace) oscillations of 
the fields around the soliton, 



hm 



^ ■ STr*id 
2 



AE TS = — ■ (Tr*K% - Tr* {K°Y 
and the second one accounts for vacuum fluctuations of the fields 

AE 

Therefore, 



^•(Tr4-Tr(A^ 



STrA| 



A m£ s 



AE- 



TS 



AEb = ^fS'I.V A- 



STr K§ 



(12) 
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Before proceeding, a further explanation of how we choose small fluctuations belonging to 
L 2 (IR 2 ) deserves a pause. Following the conventional QFT approach, we put the system in a very 
large but finite two-dimensional box and impose periodic boundary conditions on 5£,(x , x) and 
Sx{x). In Appendices II, III, and IV of Reference [22] it is shown that taking the infinite area 
limit at the end leaves no remnants for kinks and self-dual vortices. Thus, this procedure uses 
invisible boundary conditions in such a way that, at intermediate stages, one works on circles 
or genus one Riemann Surfaces. Nevertheless, the rapid (exponential) decay of the Higgs field 
to its vacuum value in these cases ensures that, starting from large lengths or areas, there will 
be a very small dependence of the results on the size. Except for the embedded ANO vortices, 
all the Higgs fields of the semi-local self-dual topological solitons decay to their vacuum values 
as for some positive k. Therefore, we expect a more significant dependence on the size of 
the box for these less concentrated solitons. 

We might try topological boundary conditions like those used in Reference [14J for the 
bosonic and supersymmetric kink. In this planar gauge theoretical setting, the analogous form 
of the anti-periodic boundary conditions would be: 

x 2 ) = e l f p ™ L ^^'^fc 1 . 5$(xi + mL, x 2 + mL) 

The line integral is along a path starting at the point P = (^1,^2), ending at P m L = (21 + 
mL,X2 + mL), and passing through regions far away from the vortex core. These twisted 
boundary conditions are also invisible in the sense that no boundary is introduced. Rather, the 
fact that we are dealing with a non-trivial line bundle of first Chern class equal to one over a 
genus one Riemann surface is taken into account. In a purely bosonic framework, however, the 
twisted boundary conditions, like periodic boundary conditions, will not leave any mark at the 
infinite area limit. 



Counter-term energies 

The Casimir energy of the previous Section is of order h, and this is also the order of the 
counter-terms found in Section 2. Thus, the one-loop semi-local topological soliton mass also 
receives contributions from the counter-terms for scalar and vector fields. At the self-dual limit 
K 2 — 1, these contributions are 

AMf, = ^[4/(l)+/(0)] J dVHSxMSal 2 ) , AM& = -hm[I(l)+I(0)] J d 2 xV k V k . 

We now reshuffle the sum of these two quantities into two pieces, respectively proportional to 
1(1) and 7(0): 

AM^ = ^/(l)£W(S,V fe ) , Ail^ = ^I(0)^(S,V k ) , (13) 

where 

^\S,V k )=4 J dPxil-^-^-^V,) , E®(S,V k ) = J d 2 x (l-\S 1 \ 2 -\S 2 \ 2 -2V k V k ) . 

The total contribution to the one-loop mass shifts from the mass renormalization counter-terms 
is: 

A M c . t . = AM c I( t 1} + AM'f 
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It follows from the classical degeneracy of vortices that there are some static deformations of a 
topological solution that do not cost energy at the tree level: those that give rise to a new self- 
dual soliton. Therefore, the dimension of the kernel of the operator K ruling the bosonic small 
fluctuations around the topological soliton should be equal to the dimension of the moduli space 
of self-dual semi-local soliton solutions. To check this, we proceed as in the previous Section by 
expanding the bosonic fields around the solution, but this time we consider only static small 
fluctuations: 

= S(x) + 5S(x) , Aj(x) = Vj(x) + 5 aj (x) 

The modified fields are still solutions of the first-order equations if (TBI) and (|T3]) are satisfied: 

F 13 = hl-^)^-d 2 6a 1 + d 1 Sa 2 + hs^SS+SS^S) = , (14) 

and 

D 1 $ + iD 2 $ = 0o(y 1 +iV 2 )5S-i(5a 1 +i5a 2 )S = , (15) 

where = dj, — iVj,. Also as in the previous Section, in order to avoid pure gauge fluctuations 
we set the static background gauge: 

+ i ( ste-^s) = o . 

Thus, the tangent space to the moduli space of self-dual vortices is the kernel of the first- 
order deformation operator D: 



-d 2 




si 


si 


si 


si 




-o 2 


-S 2 


si 


-si 


si 


s\ 


-Sf 


-d 2 + Vt 


-<9i - V 2 








si 


si 


d 1 + V 2 


-d 2 + Vi 








si 


-si 








-d 2 + Vi 


- V 2 


si 


si 








dx + V 2 


-d 2 + Vi 



One easily checks that K = V^V, and it is also possible to prove that the Kernel of T>^ is empty 
[24J. Thus, K and T> have the same kernel and the dimension of that kernel is the dimension 
of the classical moduli space: 4/. 



Zeta-function regular izat ion 

The zero-point energies for the topological soliton and the vacuum are formally super-traces: 
the differences of traces of differential operators. Such traces are divergent quantities. The 
standard procedure for dealing with this delicate point is understanding these traces as the 
generalized zeta functions of the corresponding differential operators; see [12]- [6]- [9]. Given a 
differential operator A acting on an L 2 space of functions, the corresponding generalized zeta 
function is: 

Us) = E V > see , 

Spec A 

where X n are the eigenvalues of A. We shall regularize the zero-point energies in the form: 
AE TS (s) = ^ (J^J {( K (s) - ( K c(s)} , AEo(a) = ^ {Ck (s) - Cko(s)} , 
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continuation at the point s of the s -complex plane. The physical limit is: 
AM$ S = lim AM T c s (s) = lim (AE TS (s) - AE (s)) . 



2 



Because Kq and Kq are free Schrodinger operators, their zeta functions are well known pj], [15] : 
m 2 L 2 r[s - 1] m 2 L 2 1 m 2 L 2 T[s - 1] 

Cko(s)-— f(^ + ^r' (s-i)r(s) ' flTT ' 

where T(s) and 7[s, a] are the complete and incomplete Euler Gamma functions, respectively. 

The contribution from the mass counter-terms also involves divergent quantities propor- 
tional to the integrals 1(1) and 1(0). To regularize these integrals, we apply the residue theorem 
to integrate 1(c) in the complex fc -plane. On a square of area m 2 L 2 each integral becomes a 
infinite sum over discrete momenta: 

1 f d 2 k 



_ 1 r d'k 1 _ 1 1 



2 J {I*) 2 y/k-k + c 2m2L2 £^ 2 Vk-k + c 

Accordingly, these integrals are the generalized zeta functions of the Euclidean Klein-Gordon 
and Laplacian operators evaluated at s = |. In this way, we also define the mass renormalization 
correction as a meromorphic function in the complex s-plane: 

AM«,( S ) = JL_ (C_ A+1 (s)^(S,V k ) + C-A(s)^ 0) (S,V k )) , 



and take the physical limit: 



AMft s = lim AM^ s (s) 



While — A +1 is exactly Kq, and its zeta function has been written above, calculation of Ca( s ) 
is a bit tricky. From the partition function for the Laplacian, via the Mellin transform, we 
have: 

C - A (s) = ^Urn^ (jf Wr*e-« + £ Wpr*e-*) , 
C-a(s) = ^ Km ^z^y (j[s -l,e) + T[s - 1, e)) . 



or, 



Because j[s — l,e] = ^—j- — — and T[s — l,e] = T[s — 1] , we obtain: 

, , , m2 L 2 1 

(-A « = — : t 1UV : if Re s < 1 

An (s — i)r(s) 

We finally find: 



2w?L^'^ Ly 2 J 8tt r(i) 4tt w 2m 2 L 2 ^" v 2 y 

Contrary to the kink cases, which are one-dimensional problems, a finite answer is obtained 
in the regularized integrals via the associated zeta functions. The reason is that in this two- 
dimensional problem the physical limit s = \ is not a pole of the zeta functions and only finite 
renormalizations will be necessary. 
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and generalized zeta function 



Control of Ck{s) and Ck°{ s ) is much more difficult. A convenient way for dealing with the 
zeta functions of differential operators acting on infinite dimensional spaces is by means of 
heat kernel techniques. In this Section, we shall develop this method, applied to our soliton 
operators. 



The heat kernel of a differential operator. Seeley densities 

The heat equation kernel of a N x N matrix differential operator of the general form 

K = K + Q k (x)d k + V(x) 
is the solution of the /C-heat equation 

d 



with initial condition: K)c(x, y; 0) = I-8^ 2 '(x — y). We are particularly interested in the diagonal 
x = y heat-kernel, because the Mellin transform of the partition function 



Tr e _/3/c = tr / d 2 x K K (x , x ; (3) 



gives the generalized zeta function. 
To find the kernel, one writes [8 J 

K K (x,y;(3) = C K (x,y;(3)K Ko (x,y;P) 

where /Co is the operator K, very distant from the origin, where Q k (x) and V(x) take their 
asymptotic constant vacuum values. C/c(x,y; /3) satisfies the N x iV-matrix transfer equation 

Jjl + ^j^W - ^Q k ) -A1 + Q k d k + F | C K {x , y; /3) = 

and is the unit matrix C;c{x,y; 0) = I at infinite temperature. 

Solving the transfer equation as an inverse-temperature power series expansion, 

oo 

c K (x, y-,(3) = °n(x, y; £)/? n , 

n=0 

the PDE equation becomes tantamount to the recurrence relation between the densities c n (x, y; /C): 
[nl+ (x k - y k ){d k I - ^Q k )]c n {x,y;IC) = [AI- Q k d k - V]c n „ x {x, y;JC) , n>l 

to be started from: co(x,y]K,) = I. While it is easy to find the first diagonal density, 
c\(x, x; /C) = —V(x), the determination of higher-order densities becomes more and more in- 
volved. To make the problem more tractable, we introduce the following notation: 

^C^(x) = lim d0l+a 2f d B x f ^ ^ , [c n ]M^) = {m C* B {x) . 
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densities can be written in compact form: 



(k + a x + a 2 + l)^' a ^C£* x (x) = ^ +2 ^C^ B {x) + ^ a ^C£ B {x) 

N «i 02 

-EEE 



N oil a 2 /_ \ /_ \ rqr+t-oAD 



Oil \ Oi 2 

r [t 



(ai-r+l,a 2 -t)^fDB/^\ + 



dx"[dx2 



d=l r=0 t=0 
f)r+tf)AD 

TV ai-1 a 2 / -\ / \ p.r^tr^AD 



A e e x> ( Q v (?) — >o*> + 

d=l r=0 t=0 v / \ / 12 



d=l r=0 t=0 

A 1 " «2 CKl 



zl r=0 t=0 

to be solved starting from 



w— 1 r=n ^=n \ / \ / i 2 



Cq(x, x; /C) — I =>- < (o,o)£< 



[a,fi C AB^ = o, if a ^ 0, and/or /3 ^ 



MA, 



£) = 1, A = l,2,---,iV 



The Mellin transform of the asymptotic expansion 

We must now deal with the cases N = 6 and N — 1, respectively for the operators and _fT G . 
A good approximation to the generalized zeta functions of both operators is obtained from the 
Mellin transform [§] 

1 POO 1 poo 

( K (s) = — dp^Tre-f"* , ( K c(s) = — d(3 (3 s ' 1 Tr e^ K ° 

r ( s ) Jo r ( s ) Jo 

applied to the high-temperature expansion of the partition functions 

1 oo / 4 6 

^ n=0 \ A=l A=5 



= i^E^ 



where 

W n {K)]AA = j d 2 xe-P[c n } AA (x,x;K) , A= 1,2,3,4 

[c ^)]^ = [ d 2 x[c n } AA (x,x;K) , A = 5,6 , [c„(K G )] = / rf 2 x [c n ](f , f; K G ) . 



The factor e _/3 , which appears in front of [c^if)]^ for A = 1,2,3,4, obeys the fact that the 
corresponding modes in £ have one unit of mass, while the modes for A = 5, 6 are massless. The 
generalized zeta functions are thus divided as sums of meromorphic -high-temperature regime- 
and entire -to be neglected, low temperature regime- functions of s: 
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Ck(s) 



L W n=0 [a=1 J ° .4=5 ^ J 



+ 



4vrr 
1 



d/3Tr*e 



-I3K 



+ n-l) \ T(s) 



B K (s) 



I _ °° _ /*1 -a /*oo 

Ck*00 = w^E/ d[3 f3 s+n - 2 c n (K G )e-P + — d(3Tr*e^ KC 



7[s + n-l,l] 1 



n=0 



47rrfs1 



where 7(5 + n — 1, 1] are incomplete Euler Gamma functions. We shall neglect the entire parts 
B(K) and B(K G ) and keep a finite number of terms, N , in future use of these generalized 
zeta functions for the regularization of ultraviolet divergences. 

The high-temperature one-loop semi-local vortex mass shift formula 

The contribution of the c\ coefficients to the semi-local topological soliton Casimir energy is 



AM 



(i)c. 



TS 



h ( ji 



.2 \ s 



J>i(if)U- Cl (tf ( 



A=\ 



7[M] 
47iT(s) 



but the first Seeley coefficients due to bosonic and ghost fluctuations, respectively, give: 



A=l 

X>°(#)W 



A=5 



ci(K 



d 2 x [5 - 5 1 1 2 - 3|5 2 | 2 - 2V k V k 



d 2 xll-\S 1 \ 2 -3\S 2 \ 2 -2V k V k ] 



d 2 x [1- |Si| 2 - \S 2 



Therefore, 



AM« c (-l/2) 



7[-l/2,l] 
/2) 



(^.(5,^-2/^1^).^} 



On the other hand, the contribution to the one-loop semi-local string tension shift of the mass 
renormalization counter-terms id3 

AM* (1/2) = ^ • U®{S,V h ) ■ - ■ Z<®(S } V k ) 



16tt 



r(i/2) r(i/2) 



2 7[— 1/2, 1] instead of r[— 1/2] is used in this formula to be consistent with the approximation in AM^g C (— i) 
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c\ coefficients. We finally obtain the high-temperature one-loop semi-local topological soliton 
mass shift formula: 



. , , Km 
AM TS = 



167T\/7r 



'N ( 4 6 r 

E VDfnWW - c n {K G )} ■ 7 [n - - , 1] + E 1 

n=2 I A=l A=5 



Cn(^)] 



^ ci 2 x |5 2 | 2 (xi, x 2 ) ■ f t[-^i !] 




87Tv^ J 1 1 v ' ' V 2 

In this final formula ffToT) there are four types of terms: 

• First, polynomial expressions in incomplete Gamma functions times the heat-kernel ex- 
pansion coefficients for K - saving only the four first diagonal contributions due to massive 
bosonic particles- and K G -coming from fermionic massive particles. All of them start 
from the second-order coefficients. 

• Second, polynomial expressions in — ^y, including the last two diagonal heat kernel coef- 

ficients that collect the contribution of massless Goldstone particles. The starting coeffi- 
cients are also of second order. 

• Third, a factor proportional to 41, taking into account the subtraction of the 4/ zero 
modes. 

• An extra piece proportional to the norm of the second Higgs field due to the imperfect 
cancelation of the contribution of first-order Seeley coefficients by mass renormalization 
counter-terms. 

Finally, let us mention that by cutting the expansion at a finite number, Nq, we admit an 
error - besides the rejected entire parts - which is a priori proportional to 7[-AT — ~, 1] ~ N _ \ , for 

iVo large. Nevertheless, as we shall see in the next Section, once all the calculations have been 
done the degree of convergence of the results is very good, meaning that the proportionality 
coefficient in that error is very small. The reliability of the method is therefore quite high. 



6 Mathematica calculations 

Seeley densities for spherically symmetric semi-local vortices 

We shall apply ffTBT) to spherically symmetric vortices. The heat kernel local coefficients, how- 
ever, depend on successive derivatives of the solution. This dependence can increase the error 
in the estimation of these local coefficients because we are handling an interpolating polynomial 
as the numerically generated solution, and the successive derivation with respect to r of such a 
polynomial introduces inaccuracies. Indeed, this operation is plugged into the algorithm that 
generates the local coefficients in order to speed up this process. It is thus of crucial importance 
to use the first-order differential equations Q-©-© in order to eliminate the derivatives of 
the solution and to write the local coefficients as expressions that depend only on the fields. 
Recalling the form of the spherically symmetric solutions, 

Sl(x 1 ,x 2 ) = f(r) cosO Sl(x 1 ,x 2 ) = f(r) sin (9 

S\(xi,x 2 ) = h{r) Sl{x 1 ,x 2 ) = Q 

Vi(xi,x 2 ) = sin# V 2 (xi,x 2 ) = cos# , 
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dS\ 
dx 2 
dV x 
dx\ 
dV 2 



^ [1 - a(r) cos 2 e] M = h _ a(r) sin 2 1 

dS\ _ f(r) 



f(r) 



a(r) cos 9 sin 



r 

cos 20 



-a(r)) sin 9 cos 



2/(r)«(r) / 2 (r ) + h 2 (r) - 1 



gVi 

dx 2 

„„a(r) cos 2 #.„ 9 . . , 9 . . „ . dVo 
■ 008 20^ —(f(r) + h 2 (r) - 1) — ± 



r 2 

.2 , 



cos 19 



M = _ 

dx 2 

M = _ 

ctei 
a(r) sin 



/i(r) 
r 



a(r) cos 9 
a(r)) sin 



-(f\r) + h 2 (r)-l) 



cos 20 



2 w v ' ' " v ' ' dx 2 2 
Bearing this in mind, we solve the recurrence relations to find: 



+ 



2f(r)a(r) f 2 {r) + h 2 {r) - 1 " 



tr c[(r) 
trc!,(r) 



tr4(r) 



2a(r) 2 



5f(r) 2 -3h(r) 2 



12 r 



120 r 1 



^ [4 a(r) 4 + 27 r 4 /(r) 4 - 8 r 2 a(r) (-1 + 14 /(r) 2 + fe(r) 2 ) + 
+8 a(r) 2 (-2 - 3 r 2 + 9 r 2 /(r) 2 + 3 r 2 /i(r) 2 ) + 

+ /(r) 2 (56 r 2 - 64 r 4 + 34 r 4 h(r) 2 ) + r 4 (37 - 32 h{r) 2 + 7 /i(r) 4 )] 

|-4«(r) 6 -4r 2 «(r) 3 (l4 + 35/(r) 2 - 36 h(r) 2 ) + 
r b I V / 

+4 a(r) 4 (20 + 9 r 2 + 32 r 2 /(r) 2 + 26 r 2 /i(r) 2 ) - 

-2r 2 a(r) [57 r 2 /(r) 4 + /(r) 2 (32 + 331 r 2 - 75 r 2 /i(r) 2 ) - 4 (-1 + h(r) 2 ) (-16 - 9r 2 + r 2 h(r) 2 )] + 
+a(r) 2 [-256 - 144 r 2 - 117 r 4 + 99 r 4 /(r) 4 - 16 r 2 h(r) 2 + 94 r 4 h(rf - 61 r 4 h(r) 4 + 

+ 2r 2 /(r) 2 (56 + 183 r 2 + 19 r 2 /i(r) 2 )] + 
+r 4 [-16 + 151 r 2 - 29 r 2 f(r) 6 + (32 - 135 r 2 ) h(rf + (-16 + 23r 2 ) fe(r) 4 + r 2 h(rf + 

+ /(r) 4 (-20 + 199r 2 -57r 2 /i(r) 2 ) + /(r) 2 (392 - 321r 2 + 2 (-68 + 111 r 2 ) h(r) 2 - 27 r 2 /i(r) 4 )] } 



trcf(r) = i_^JL_/( r )a_ 3ft ( r )3 

tr^(r) = -i-[4a(r) 4 -r 4 /(r) 4 + 8r 2 a(r) (l + 2/(r) 2 -A(r) 2 ) - 

-8a(r) 2 (2 + r 2 + r 2 /(r) 2 -5r 2 fe(r) 2 ) +2r 2 /(r) 2 (-4 + 9 r 2 A(r) 2 ) + 
+r 4 (l - 8 h(rf + 19 ft(r) 4 )] 

trc?(r) = ^ {4a(r)« - 4r 2 a(r) 3 (-14 + 9 /(r) 2 + 84Mr) 2 ) - 
-4a(r) 4 (20 + 3r 2 + 2r 2 (/(r) 2 + 4 /i(r) 2 )) + 

+a(r) 2 [256 + 48 r 2 -3 r 4 + 45 r 4 /(r) 4 + 2 r 2 (-40 + 89r 2 ) h(r) 2 - 115 r 4 h(r) 4 + 

+ 2r 2 /(r) 2 (8 + 5r 2 - 35 r 2 /i(r) 2 )] - 
-2r 2 a(r) [53 r 2 /(r) 4 + 4 (-1 + /i(r) 2 ) (-16 - 3 r 2 + 7r 2 ft(r) 2 ) - /(r) 2 (32 + 17r 2 + 47 r 2 fe(r) 2 )] + 
+r 4 [l6 + 3r 2 +3r 2 /(r) 6 - (32 + 19r 2 ) h(r) 2 + (16 + 23 r 2 ) h(r) 4 + 33 r 2 h(rf + 

+ /(r) 4 (52-r 2 + 39r 2 /i(r) 2 ) + f(r) 2 (-24 - 5 r 2 + (-72 + 22 r 2 ) /i(r) 2 + 69 r 2 h(r) 4 )] } 



cf(r) = l-/(r) 2 -Mr) 2 

c G( r ) = _L [3r 2 +2r 2 /(r) 4 - (5r 2 + 4«(r) 2 | h(rf + 2 r 2 h(r) 4 + 
+fif) 2 (- 4 - 5r 2 +8a(r) -4a(r) 2 +4r 2 M r ) 2 )] 
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+r 2 [-17r 2 + 8a(r) - 16a(r) 2 ] /i(r) 4 +4r 4 fc(r) 6 + 

+r 2 /(r) 4 [-24 - 17 r 2 + 40 a(r) - 16 a(r) 2 + 12 r 2 /i(r) 2 ] + 

+/( r ) 2 [-32«(r) 3 + 16a(r) 4 + 8r 2 a{r) (-5 + 6 h(rf^j + 16 a(r) 2 (l + r 2 - 2r 2 /i(r) 2 j + 
+ r 2 ^24 + 23r 2 - 2 (10 + 17r 2 ) /i(r) 2 + 12 r 2 M' r ) 4 )] } 

etcetera, by means of a computing program implemented on a PC with Mathematica. 

One-loop mass shift for the mass of the ANO vortex 

Denoting simply 



trc I n = 2n / drrtrc^r) 
Jo 



trc 



o 



2rr 



dr rtrc^(r) 



2rr 



drrc^(r) 



by plugging the NO vortex solution in these expressions, i.e. the case ho = 0.0, embedded in 
this model, we find the Table at the left: 





ho = 0.0 


n 


trC n 


tr# 




1 


-41.4469 


-91.8429 


12.599 


2 


30.3736 


0.96286 


2.61518 


3 


12.9447 


-0.0592415 


0.32005 


4 


4.22603 


0.001512548 


0.0230445 


5 


1.05059 


0.000758663 


0.0013023 


6 


0.20900 


-0.00023912 


0.0000698185 



N 


AM v (No) 




I = 1 


2 


-1.61536 


3 


-1.66862 


4 


-1.67809 


5 


-1.67966 


6 


-1.67989 



(left) The sixth lowest Seeley coefficients for Nielsen- Oles en self-dual vortices, (right) 
Convergence of the one-loop mass shift for self-dual semi-local NO vortices in units of h~m. 



Because the NO vortex solutions have been generated numerically, integration over the whole 
plane of the Seeley densities can also only be performed numerically. Therefore, we are forced 
to put a cut-off into the area and replace the infinite plane by a discus of radius R, which in the 
calculations above was chosen to be R = 10.000 (compare with the profiles of Figure 1). Use 
of these numbers in formula (f!6l) provides us with the Table at the right, where the one- loop 
mass shifts in km units of semi-local self-dual NO vortices are shown up to sixth order in the 
asymptotic formula. 

Our result for the one-loop mass shift of semi-local self-dual NO vortices is: 



AM V (N = 6) = -1.67989ftm 



(17) 



The ratio between the mass shifts of self-dual NO vortices in the semi-local and normal Abelian 
Higgs model is: 

AM-hm _ L67989 _ 



AM™ 1.09449 

see [15J. Similar relations exist between ratios of kink mass shifts in the A(0)^ model and the 
BNRT model; a two-field system that depends on a real positive parameter a: 



am bnrtm ^ o, 6 93943 
AMf M ~ 0.471113 



1.47299 if a = 0.99 
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am a*m 0.471113 

see [IT]. 

Several comments are in order: 

1. The coefficients trc^, for n = 2, 3, • • • , 6, are identical, within numerical precision, to the 
same coefficients for self-dual NO vortices in the Abelian Higgs model. Therefore, the 
difference of the vortex mass shift arising in the semi-local Abelian Higgs model is due to 
the contribution of the tr c° coefficients and to the double number of zero modes. 

2. trc{ and trdp are very large negative numbers. In the case at hand, h Q = 0, their 
contribution cancels against the energy induced by mass renormalization counter-terms. 
A similar behavior of ho ^ self-dual semi-local topological solitons would mean that 
rather than decreasing their mass should grow as the result of one-loop fluctuations, a 
possibility that we shall study below. 

Seeley coefficients for semi-local topological solitons: 
one quantum of magnetic flux and various values of Hq. 

To accomplish this goal, we offer several Tables in the Appendix with similar numerical calcu- 
lations on discuses of radius R = 10 2 , R = 10 3 , R = 10 4 , R = 10 5 , and R = 10 6 for the values 
ho = 0, ho = 0.1, ho = 0.3, h = 0.6, and ho = 0.9. Inspection of these Tables raises some 
points that merit comment: 

1. We see from the Tables for the ho = 0.0 case that although the first-order coefficients grow 
spectacularly with the radius, the highest coefficients are quite stable against growing 
areas of integration, suggesting good behavior -a finite value very close to the value 
obtained at R = 10 4 - of the one-loop correction in the infinite area limit. 

2. Things start to be different when we consider the ho = 0.1 case. First, the positive and 
negative contributions of the first-order Seeley coefficients are not completely canceled out 
by mass renormalization counter-terms. Second, of special importance is the departure 
of the values of trc°(K) from the same numbers for the NO vortex, whereas trc^(i^) 
and c n (K G ) change slowly with h . Moreover, trcf (K) for h = 0.1, besides not being 
exactly canceled, is negative and large. This effect prompts a increasingly less negative 
one-loop correction, eventually changing sign, with larger and larger R. This behavior 
differs completely from that of embedded NO vortices and suggests a different fate for 
semi-local topological solitons when quantum fluctuations enter the game. 

3. The numbers offered in the Tables for ho = 0.3 follow a similar pattern to those of the 
h = 0.1 semi-local topological soliton. There are, however, some quantitative differences. 
trc^(i^), n ^ 1, Hq — 0.3, differs from the same number of the NO vortex more than the 
h = 0.1 topological soliton. Thus, these numbers grow with ho and R. trc[(K) decreases 
with respect to the same number for the NO vortex and with R. There is a change of sign 
at R = 10 3 if ho = 0.1 and at R = 10 4 if ho = 0.3. The numbers trcf (K) behave as in the 
h = 0.1 case but depart more from NO values. If h = 0.1 the ghost coefficients c n (K G ) 
are slightly larger than these numbers for h = 0, but they start to decrease (rapidly with 
R) when h = 0.3. 

4. The pattern observed for the evolution of the Seeley coefficients when the area grows in 
the cases h = 0.1 and h = 0.3 is reinforced for broader semilocal topological defects 
with h = 0.6 and h = 0.9. The unbalanced first-order coefficients trcf, obtained by by 
mass renormalizations, rapidly tend toward huge negative values. 
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The consequence is as follows: whereas the one-loop mass shift of embedded ANO vortices 
is always negative and varies extremely slowly as the area increases towards more negative 
values, one- loop mass shifts of genuine semilocal topological solitons with \h \ > become less 
negative, and even positive, for larger areas, as is shown in the following Table. 



One-loop mass shifts for semi-local topological solitons: Five values of ho, five values of R, and fixed N = 6. 



R 


AM V (N = 6,R) 


AM V (N = 6,R) 


AM y (iV = 6,R) 


AM V {N = 6,R) 


AM V {N = 6,R) 




h = 0.0 


h = 0.1 


h = 0.3 


h = 0.6 


h = 0.9 


10' 2 


-1.67955 


-1.61672 


-1.05000 


2.10142 


24.6066 


10 3 


-1.67971 


-1.58311 


-0.626167 


4.5485 


42.7747 


10 4 


-1.67989 


-1.55133 


-0.252586 


6.41655 


60.9433 


10 5 


-1.68005 


-1.51957 


0.12086 


8.5741 


79.1116 


10 6 


-1.68026 


-1.48779 


0.49433 


10.7203 


97.2798 



The classical degeneracy in energy between semi-local topological defects seems to be bro- 
ken by one-loop fluctuations, the embedded ANO vortices becoming the ground states in the 
topological sector of one quantum of magnetic flux. It is remarkable how strong this effect 
becomes for large \h \. 

These numerical results find support in the following qualitative arguments based in the 
analysis of the potentials arising in the matrix Schrodinger operators governing small fluctua- 
tions around semi-local self-dual topological solitons. The more pertinent, diagonal, operators 
in K and K G are: 

A = K G = - A +f 2 (r) + h 2 (r) = - A +A(r) 

1 n 2 (r) 
B = - A +-(3/ 2 (r) + h\r) + - 1) = - A +B(r) 

1 n 2 (r) 
C — — A +-(/ 2 (r) + 3h 2 (r) + 2—^ - 1) = - A +C(r) . 

A glance at Figure 2 reveals: 

1. For ANO vortices, h = 0.0 and h(r) = 0, there are attractive potential wells with the 
bottom at the origin. Thus, the correction must be negative. Also, the field profiles reach 
their vacuum values exponentially and the corrections are quite stable with respect to the 
area of the normalization box. 

2. For h = 0.1, h(r) is not zero, see also Figure 1 and take into account that: ~ r ->o r 2 . 
We observe that the corresponding terms push the wells upwards from the bottom, making 
them less attractive. This explains why the mass shifts are less negative. 

3. For ho = 0.3, h(r) is big enough to globally produce a change from attractive to repulsive 
potential forces provided that the area of the normalization box is sufficiently large. There 
is consequently a change in the sign of the mass correction. 

4. For ho = 0.6 and ho = 0.9, potential barriers dominate starting from relatively small 
areas of the normalization box. The field profiles go to their vacuum values very slowly 
and the barriers become very wide, explaining the strong dependence on the size of the 
normalization box. 

We conclude that the classical degeneracy is broken by one-loop fluctuations. Even if the 
classical topological bound is saturated by all the solitons in the moduli space of solutions of 
the first-order equations, quantum effects can distinguish between the different energy densities 
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a) Functions A{r), B(r) and C(r) for a) ho = 0.0, a) ho = 0.3 and a) ho = 0.9 



of these extended objects. The two extremes are the ANO vortices, ho = 0, where the energy 
density is concentrated around the zero of the Higgs field, and the CP Mumps, with energy 
densities uniformly distributed over the whole spatial plane. A similar effect has been observed 
before in the moduli space of degenerated two-component kinks analyzed in Reference |17j . 



Infrared divergences: quantum fate of semi-local topological solitons 

The dependence on the area of the normalization box is due to the slow decay (non-exponential) 
to their vacuum values of genuine semi-local topological solitons as compared with ANO vor- 
tices. Plugging the asymptotic form of the spherically symmetric topological soliton solutions 
obtained in Section §. 1.7 into the Seeley densities, we find the following behavior at infinity 
in terms of the parameter l/i 1 ! (which sets the large r behavior of the solutions): 

2vrrtrc((r) ~ - — (1 - \h x \ 2 ) + -^(^/i 1 ) 2 - (/i 1 ) 4 ) + 0(4) 

27rrtrcf(r) ~ - — (1 + \h l \ 2 ) + %{A\h l \ 2 + \h l \ A ) + 0(4) 

ry y OO T* <y>J rr»%J 



8tt,,i,9 1 

oo r 3 



27rrcf(r) ~ ^|/ l 1 | 2 + 0(^) 



2-tr4(r) * r ^^^\hr + ^(-l + ^\ 2 -\hr) + 0^) 
2*rtr<£(r) * f ^ ^ ^| 2 + J(-l + ^\ 2 - \h^) + 0(1) 
2*rc£(r) - r ^ oc ^ 1 | 2 + ^(7681^-801^) + ^) 

2-trc|(r) - r ^ 00 |l^l 2 + ^ 1 | 2 -|^l 4 ) + ^) 

2*rt«£(r) ^ r ^ ^ - ^r\h X \ 2 + ^{-^\h l \ 2 + (/i 1 ) 4 ) + O(^) 

2vrrc 3 G (r) ~ ^l^ + O^) 

t — > oo 5r r 8 

2rtci(r) * ^l^ 1 ! 2 + ^(^Vl 2 - l^l 4 ) + 0(4) 

r — > oo or or 13 5 r° 

2rtc?(r) ~ -^1^12 + JL(4|^|2- |^|4) + o(l) 

r — > oo or or 13 r° 

2irrc?(r) ~ ^^Ife 1 ! 4 - — l/t 1 ! 6 + — l/t 1 ! 8 + - I^P + Of— ) 



2 — - ^. . oi |ftl|2+ 3P (2|/il|2 - |/ll|4)+ ^ ) 



r 



2vrrtrc?(r) ~ - —\h l \ 2 + -^{-A\h l \ 2 + \h l \ A ) + 0(4) 

51 ; r^oo 30r' 1 30r 3 V 1 11 1 ; K r 5 ' 



r, s /2048 llll2 416,, 1|4 43 ll1lfi 1 .,,.„ 1 ,,i,irA 1 . 

2vrrc?(r ~ - /i 1 2 /i H h 1 6 /i + /i 1 10 + 0(— r- 

5W r^oo r 9 V 105 1 1 189 1 1 540 1 1 864 1 1 138240 1 1 J V 11 ' 
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2vrrtr C6 / (r) ~ *\ h i\* + *(l\ h i\2-h h i\4) + (±r 

by ' r^oo 180r' 1 15r 3 V 1 12 1 1 ; V 5 

2vrrtrc 6 °(r) ~ *\ h if + *Q h if -htff) +0 (\) 

r -> oo 180r 45r 3 4 r 5 



The key observation is the appearance of infrared logarithmic divergences in the Seeley coef- 
ficients trc^(K) and trc°(K) for all n. The ghost coefficients c^(K G ), however, are infrared 
convergent. The combination of the signs that we have seen in the previous sub-Sections and 
the large r behavior show that one-loop mass shifts of semi-local topological solitons tend to 
+oo at the infinite area limit. Semi-local topological defects grow infinitely massive due to the 
infrared effects of one-loop fluctuations. This phenomenon seems to be amazingly close to the 
non-existence of Goldstone bosons in (l+l)-dimensions. 

There is a very important exception: for ANO vortices, {h 1 ] = and only the first-order 
coefficients are infrared divergent. However, the contribution of these coefficients is totally can- 
celed by mass renormalization counter-terms. Our results suggest that only the ANO vortices 
between all the semi-local topological solitons survive one-loop quantum fluctuations. It would 
be very interesting to try a more analytic approach to this problem in order to fully elucidate 
this delicate issue. 



Appendix 



Semilocal strings with h = 0.0. Embedded NO vortices 

Seeley coefficients: Semilocal strings with ho 



0.0 





h = 0.0 , R = 10 2 




h = 0.0 , R= 10 3 


n 


tr c n 


tr c G 


c° 


n 


tr4 


trc^ 


°n 


1 


16.3087 


-33.9956 


12.5761 


1 


-12.5691 


-62.9193 


12.5875 


2 


30.3548 


0.959353 


2.61136 


2 


30.3641 


0.960952 


2.61327 


3 


12.9428 


-0.0596236 


0.319668 


3 


12.9437 


-0.0594326 


0.319859 


4 


4.2259 


0.00149819 


0.0230172 


4 


4.22596 


0.00151183 


0.0230309 


5 


1.050558 


0.000757146 


0.00122871 


5 


1.05058 


0.0000747904 


0.00122947 


6 


0.209003 


-0.000239188 


0.0000697495 


6 


0.209003 


-0.000239154 


0.000069784 




h = 0.0 , R= 10 4 




h = 0.0 , R = 10 5 


n 


tr< 


tvc° 




n 


tr< 


trc° 




1 


-41.4469 


-91.8429 


12.599 


1 


-70.3247 


-120.767 


12.6105 


2 


30.3736 


0.96286 


2.61518 


2 


30.3829 


0.964771 


2.61709 


3 


12.9447 


-0.0592415 


0.32005 


3 


12.9456 


-0.0590504 


0.320242 


4 


4.22603 


0.001512548 


0.0230445 


4 


4.2261 


0.00153913 


0.0210582 


5 


1.05059 


0.000758663 


0.0013023 


5 


1.05059 


0.000759469 


0.00123099 


6 


0.20900 


-0.00023912 


0.0000698185 


6 


0.209003 


-0.000239086 


0.00002698506 





h = 0.0 , R= 10 b 


n 


trc^ 


trc^ 




1 


-99.2028 


-149.690 


12.6220 


2 


30.3920 


0.9666682 


2.61901 


3 


12.9466 


-0.0588602 


0.320432 


4 


4.22617 


0.00155181 


0.02300754 


5 


1.05058 


0.0007600532 


0.00123155 


6 


0.209003 


-0.00023874 


0.0000695093 
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Seeley coefficients: Semilocal strings 





ho = 0.1 , i? = 10 2 




h = 0.1 , i? = 10 3 


n 


tr 






n 


tr c n 


tic n 




i 
i 


ly.uooy 


oc 1 Q«7 


IZ.oUOl 


i 
l 


7 fi7^q 
- 1 .D i OOO 


-D4.D4Zy 


lo.uoyo 


Z 


qi ooqo 
ol.ZooZ 


i.oyuoz 


9 R91 1 7 
Z.OZ11 / 


Z 


qi Q9i c; 

oi.yzio 


9 i /inno 
z. i4uuy 


Z.D0O44 


q 


lo.ZUUU 


n 977/1 1 7 
-U.Z I /41 1 


n qi P.1'7 


Q 
O 


i q q^77 


-U.4UDOUO 


n q9i qfi/i 


A 
4 




U.OU / Z4L) 


D 099^497 

VJ.UZZU4Z / 


A 
4 


A qi 


u.uyuoyvj 


U.UZZOOUO 


z. 
o 


1.06171 


-0.0102584 


0.00119161 


c; 
(j 


1.0685 


-0.01694523 


0.00121037 


6 


0.210813 


0.00159792 


0.0000667653 


6 


0.211935 


0.00271637 


0.000067618 




h = 0.1 , J? = 10 4 




/jo = 0.1 , R = 10 5 


n 




tr# 




fi 


trc n 


tr# 




i 
i 


-34.3878 


-94.099 


13.3734 


i 
i 


-61.1001 


-123.555 


13.6571 


9 

z 


32.5602 


2.58966 


2.71571 


9 
Z 


33.1989 


3.03924 


2.76299 


q 


13.5154 


-0.536179 


0.326091 


q 


13.6732 


-0.665552 


0.330819 


A 
4 


4.35407 


0.124453 


0.023218 


A 
4 


4.38928 


0.158316 


0.0235856 


5 


1.0753 


-0.0236316 


0.00122913 


5 


1.0821 


-0.0303191 


0.00124789 


6 


0.213057 


0.00383473 


0.0000684707 


6 


0.214679 


0.0049531 


0.0000693256 





ho = 0.1 , R = 10 b 


n 


trc^ 


tr# 




1 


-87.8124 


-153.011 


13.9407 


2 


33.8375 


3.48881 


2.81037 


3 


13.8312 


-0.784926 


0.335521 


4 


4.42449 


0.192179 


0.00238892 


5 


1.0889 


-0.030047 


0.00126624 


6 


0.215306 


0.00607261 


0.0000689434 



Semilocal strings with h = 0.3 



Seeley coefficients: Semilocal strings 





h = 0.3 , R= 10 2 




h = 0.3 , R= 10 3 


n 


trt£ 


tr# 




n 




tr# 




1 


27.8909 


-48.6659 


11.7815 


1 


2.4708 


-86.5124 


10.8818 


2 


39.9546 


7.69438 


2.22326 


2 


40.2110 


11.5509 


2.07325 


3 


15.2219 


-2.2614 


0.247514 


3 


16.4824 


-3.61196 


0.232514 


4 


4.77643 


0.540926 


0.0159336 


4 


5.10495 


0.873748 


0.0148622 


5 


1.15832 


-0.10716 


0.000764659 


5 


1.22480 


-0.173998 


0.000705138 


6 


0.226592 


0.0176963 


0.0000395987 


6 


0.237707 


0.0288234 


0.0000368937 




/j = 0.3 , i? = 10 4 




h = 0.3 , R= 10 5 


n 


tr4 


tr# 




n 




trc^ 




1 


-22.9489 


-124.360 


9.98213 


1 


-48.3715 


-162.207 


9.08246 


2 


43.4678 


15.4073 


1.92331 


2 


46.7241 


19.2687 


1.77336 


3 


17.7429 


-4.96241 


0.21752 


3 


19.0034 


-6.31286 


0.202525 


4 


5.43346 


1.20654 


0.0137912 


4 


5.76197 


1.53933 


0.0127201 


5 


1.29128 


-0.24083 


0.000645635 


5 


1.35776 


-0.30766 


0.000586148 


6 


0.248822 


0.0399494 


0.0000341891 


6 


0.259938 


0.0510755 


0.0000314752 



26 



n 


trc n 


tr c n 


°n 


1 


-73.7864 


-200.055 


8.18278 


2 


49.9803 


23.1201 


1.62347 


3 


20.2636 


-7.6633 


0.187535 


4 


6.09059 


1.87212 


0.011658 


5 


1.42424 


-0.374493 


0.00052486 


6 


0.271056 


0.0622017 


0.0000284767 



Semilocal strings with h = 0.6 



Seeley coefficients: Semilocal strings 





h = 0.6 , # = 10 2 




h = 0.6 , R = 10 3 


n 




tr# 


a 


n 




tr# 




1 


107.840 


-113.330 


13.3675 


1 


132.751 


-189.607 


14.4510 


2 


72.5374 


41.052 


1.73758 


2 


97.6535 


65.4471 


1.91781 


3 


26.5220 


-12.8641 


0.135265 


3 


34.6838 


-20.9178 


0.153289 


A 

4 


7.45753 


3.13219 


0.00625589 


4 


9.48185 


5.15141 


0.0075433 


5 


1.67277 


-0.615441 


0.000232254 


5 


2.0767 


-1.01896 


0.000303781 


6 


0.309679 


0.101414 


8.66613x10^ 


6 


0.376957 


0.168681 


0.0000119181 




ho = 0.6 , R = 10 4 




h = 0.6 , R = 10 5 


n 


trc n 


trc T V 


a 


n 




trc^ 




1 


157.661 


-265.885 


9.98213 


1 


182.571 


-342.162 


16.6184 


2 


122.770 


89.8409 


1.92331 


2 


147.886 


114.235 


2.27905 


3 


42.8451 


-28.9709 


0.21752 


3 


51.0065 


-37.0239 


0.189413 


4 


11.50061 


7.17046 


0.0137912 


4 


13.5303 


9.18952 


0.0101236 


5 


2.48061 


-1.42244 


0.000645635 


5 


2.88453 


-1.82592 


0.000447104 


6 


0.44232 


0.235944 


0.0000151754 


6 


0.511506 


0.303206 


0.0000184335 





h = 0.6 , R = 10 s 


n 




tr c^ 


°n 


1 


207.481 


-418.440 


17.7021 


2 


173.003 


138.629 


2.45961 


3 


59.1167 


-45.0769 


0.207444 


4 


15.5545 


11.2086 


0.0114354 


5 


3.28844 


-2.2294 


0.000517096 


6 


0.578786 


0.370467 


0.0000217957 



Semilocal strings with h = 0.9 

Seeley coefficients: Semilocal strings 
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n 


tr< 


trc- 


r 

°n 




n 


trc n 


trc- 




1 


638.887 


-601.992 


17.3563 




1 


1055.500 


-1030.530 


25.0050 


2 


326.423 


290.227 


1.24111 




2 


536.449 


495.158 


2.5130 


3 


109.814 


-94.6224 


0.0897853 




3 


178.399 


-162.381 


0.216979 


4 


27.9842 


23.3075 


0.00588282 




4 


45.0018 


40.2880 


0.014969 


5 


5.71851 


-4.59264 


0.000325247 




5 


9.11549 


-7.798641 


0.00083022 


6 


0.97486 


0.755329 


0.0000152661 




6 


1.54073 


1.32206 


0.0000282526 




h = 0.9 , R = 10 4 




ho = 0.9 , R = 10 5 


n 


tre£ 


tr# 




n 




trc^ 




1 


1472.190 


-1459.110 


32.6701 


1 


1888.870 


-1887.670 


40.3354 


2 


746.482 


700.081 


3.79044 


2 


956.515 


905.004 


5.06804 


3 


246.919 


-230.135 


0.344728 


3 


315.440 


-297.889 


0.472483 


4 


62.0179 


57.2676 


0.0240939 


4 


79.0340 


74.2472 


0.0332125 


5 


12.5121 


-11.3800 


0.00133716 


5 


15.2088 


-14.7736 


0.000184412 


6 


2.10653 


1.88677 


0.0000612957 


6 


2.67232 


2.45247 


0.0000843383 





h = 0.9 , R= 10 6 


n 




trc^ 




1 


2305.550 


-2316.230 


48.007 


2 


1166.550 


1109.930 


6.34561 


3 


383.961 


-365.644 


0.60024 


4 


96.05501 


91.2268 


0.0423323 


5 


19.3054 


-18.1672 


0.0423323 


6 


3.23313 


3.01818 


0.000107346 
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